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A semi-analytical approach to atomic waveguide scattering for harmonic confinement is developed 
taking into account all partial waves. As a consequence ^-wave confinement-induced resonances are 
formed being coupled to each other due to the confinement. The corresponding resonance condition 
is obtained analytically using the if-matrix formalism. Atomic scattering is described by transition 
diagrams which depict all relevant processes the atoms undergo during the collision. Our analytical 
results are compared to corresponding numerical data and show very good agreement. 
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I. INTRODUCTION 

Reduced dimensionality in degenerate ultracold atomic 
gases plays a key role for the experimental realization 
and theoretical investigation of exotic quantum phases 
such as the Tonks- Girardeau gas Specifically, atomic 
collisions in waveguides represent a fundamental ingredi- 
ent for studying the effects of the reduced dimensional- 
ity. Indeed, in two-body collisions, the confinement im- 
plies significant modifications of the scattering proper- 
ties of the collisional complex. Existing theoretical stud- 
ies on bosonic collisions show that resonant scattering 
can be induced by the confinement, yielding the so-called 
confinement-induced resonance (CIR) effect [2|,l3j . A CIR 
emerges when the length scale of the confinement be- 
comes comparable to the s- wave scattering length of the 
colliding bosons and it is interpreted as a Fano-Feshbach- 
like resonance. An important characteristic property is 
that the effective one-dimensional two-body interaction 
can be controlled by adjusting the confinement parame- 
ters. 

In recent years, the significant progress of quantum 
technologies has lead to the experimental observation 
and exploration of the CIR physics in quasi-one- and 
quasi- two- dimensional waveguides for both, bosons [J- 
0, [1| and fermions 0, 0] • Additionally, remarkable the- 
oretical studies focus on CIR physics and its extensions 
such as p-wave CIR in spin-polarized fermions llOL dual- 
CIR [ll| , and the resonant molecule formation [12j . Fur- 
thermore, the concept of CIR physics has been investi- 
gated also in the case of quasi- two-dimensional or an- 
harmonic waveguides jl3l4l8| and in the case of multi- 
channel or mixed-dimension scattering (191421}. Higher- 
partial wave interactions constitute an interesting exten- 
sion within the framework of CIR physics as it was shown 
in [22j, since they are expected to provide novel many- 
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body phenomena, such as unconventional superconduc- 
tivity or superfluidity [23l-[26|. Nevertheless, a thorough 
theoretical treatment is still needed extending the CIR 
physics beyond s- or p-wave interactions. 

In this work, we present a theoretical framework 
for atom-atom collisions in the presence of a harmonic 
waveguide, where all the higher partial waves are prop- 
erly taken into account either for bosons or fermions. 
Mainly, it consists of a generalization of the if-matrix 
approach for quasi-one-dimensional geometries presented 
in [lOj. Our treatment yields mutually coupled £-wave 
confinement-induced resonances as well as an analytical 
relation for the position of all ^-wave CIRs. Further- 
more, converting the "physical" quasi- ID K- matrix into 
a Dyson-like form of equations allows us to classify, in 
terms of specific transition diagrams, all possible pro- 
cesses which the two atoms undergo during the collision. 
This detailed analysis permits us to illustrate the con- 
nection of the transitions in the manifold of the closed 
channels of the transverse confinement with the Fano- 
Feshbach scenario for resonant scattering. The validity 
of our method is verified showing an excellent agreement 
of the obtained analytical results with numerical calcu- 
lations of the transmission coefficient in a system of two 
bosons interacting via a Lennard- Jones potential in the 
presence of harmonic confinement. 

Our work is organized as follows. In Sec. II we discuss 
the general aspects of K-matrix theory, and present the 
if-matrix approach for quasi-one dimensional systems. 
Sec. Ill contains a comparison of the analytical results 
with the numerical calculations and a corresponding dis- 
cussion. Sec. IV provides a brief summary and conclu- 
sions. 



II. if-MATRIX REPRESENTATION FOR 
QUASI-ONE DIMENSIONAL SYSTEMS 

The if-matrix is defined as a Caley transformation of 
the ^-matrix given by the following relation: 
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1-iK' 
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The form of eq.(pQ) is such that the hermiticity of the ma- 
trix K ensures the unitarity of S. The advantage of the 
If -matrix with respect to the ^-matrix is that it possesses 
a simpler pole structure. This allows for an analysis of 
overlapping resonances in a straightforward manner. 

Let us now consider collisions of two atoms in the pres- 
ence of an external potential, which confines the transver- 
sal degrees of freedoms. The transversal confinement is 
induced by a harmonic potential, which permits a sepa- 
ration of center of mass and relative degrees of freedom. 
Thus the two-body scattering problem can be reduced to 
a problem of a single particle (relative particle) scattered 
by a potential (interatomic potential) placed at the ori- 
gin. The Hamiltonian of the relative particle, expressed 
in cylindrical coordinates, reads: 
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pu 2 ± p 2 + V^(r), 
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where /i is the reduced mass of the two colliding atoms, 
uo±_ is the transversal confinement frequency and V^(r) 
is a short-range interatomic potential. We assume that 
the interatomic and the transversal confinement potential 
do not depend on the azimuthal angle <p. Consequently, 
the corresponding quantum number m is conserved and 
throughout this paper we will focus on m = 0. 

The relative Hamiltonian imposes two approximate 
symmetries in different regions of the configuration space. 
Near to the origin, r < ro, the interatomic potential 
Vsh(r) of range ro dominates and therefore the collision 
process obeys spherical symmetry. Thus in this regime, 
for r « ro, the relative particle experiences a free-space 
collision off the interatomic potential with the total col- 
liding energy E = h 2 k 2 /2p. This process can be de- 
scribed by the following radial wavefunction: 



(3) 



where a labels the linearly independent solutions, i^(r, 0) 
(Ge(r,0)) is the regular (irregular) solution, i.e. spheri- 
cal Bessel ji(r) (spherical von Neumann ne(r)) functions 
multiplied by the Legendre polynomials Pi{6), and the 
summation runs over all the even (odd) £ for the case of 
bosons (fermions). Kf® are the elements of the if- matrix 
K 3D in three-dimensions encapsulating all the scattering 
information related to the interatomic potential V^(r). 
Due to the short-range character of the interatomic in- 
teractions K_ 3D is diagonal and its elements are equal to 
tan 5^, with Si being the phase shift of the ^-th partial 



Additionally, in the asymptotic limit, |r| — >> oo, the 
transversal confinement term in eq.(j2j) dominates imply- 
ing its cylindrical symmetry. The total colliding energy 
E can be written as a sum of the transversal energy being 



determined by the harmonic confinement and the longi- 
tudinal energy, which involves to the z-component of the 
momentum of the relative particle, according to the re- 
lation E = H 2 k 2 /2p = Huj ± (2n + 1) + H 2 q 2 /2p. In this 
region the wavefunctions are axially symmetric and can 
be written as follows: 
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g n (z,p)K} D 
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where f3 labels the linearly independent solutions, 

are the elements of the if-matrix K_ 1D in quasi-one di- 
mension and f n (z,p) (g n (z,p)) is the regular (irregular) 
solution in the presence of the trap with no interactions. 
The specific form of the regular and irregular solutions 
are as follows: 
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for bosons 



for fermions 



cos q n \z\ 
sin q n \z\ 
_z_f sin<? n |z| 
M ^ - cos<? n |z| 

(5) 

where $> n (p) are the eigenfunctions of the two- 
dimensional harmonic oscillator for m = and the func- 
tions in the curly bracket refer to the longitudinal de- 
gree of freedom including the symmetrization (antisym- 
metrization) for the case of bosons (fermions). 

So far it is clear that the scattering process of atomic 
collisions in the presence of a harmonic waveguide is 
characterized by two distinct regions in the configura- 
tion space which possess different symmetries. Thus, the 
key idea is to perform a frame transformation which will 
permit the matching of eqs.© and (|4]) yielding an expres- 
sion of K_ 1D in terms of K^ D . Obviously, a unitary frame 
transformation in this case doesn't exist since eqs.(|3J) and 
(HJ do not fulfill the same Schrodinger equation in com- 
plete configuration space. However, the length scales of 
the inte ratomic potential ro and the harmonic oscillator 
a± = ^h/puj± differ by orders of magnitudes, ro <C a± 
This means that between the two distinct regions 
of spherical and cylindrical symmetry there is a region 
where V s h(r) w ^puo\p 2 w 0. Consequently, for this con- 
figuration subspace eqs.© and (|4]) fulfill approximately 
the same Schrodinger equation, and one can perform a 
local and non-orthogonal frame transformation. The idea 
of local frame transformations was introduced by Harmin 
and Fano [2?], [HI and extended later on by Greene [29| . 

As it was shown in [29| the above-mentioned local 
frame transformation does not depend on any additional 
parameter in order to match the two sets of solutions, 
eqs.© and (|4]). The only condition which has to be 
fulfilled is ro <C a± (separation of scales) yielding the 
following relations, which connect eqs.© and (|4]): 

f n (z, p)= ^2 ^( r > 0) u tn and 



g n (z,p)=Y,Gi(r,0)(U T y £ 



(6) 
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where Ue n is the matrix element of the local frame trans- 
formation U given by the relation: 



V2(-l) do 2£ + l 



a± 



kq n 



(7) 



In eq. (O Pi(^r) is the ^-th Legendre polynomial and the 
term do is £/2 or {£ — l)/2 for the case of bosons or 
fermions, respectively. Note that the local transforma- 
tion U connects the £ partial waves with the n-modes of 
the transversal confining potential. Thus K_ 1D with the 
help of eqs. (|4]) and © will take following form: 



K 1D = U T K 3D U. 



In eq. (j8]) the matrix K involves an admixture of all n 
channels of the transversal confinement with all £ partial 
waves of the interatomic potential for all energies. 

In our study we focus on the single-mode regime, 
where the collision takes place in the ground state of 
the transversal confinement. Thus, the total colliding 
energy is E = h 2 k 2 /2p = Huj± + T?q^j2\i yielding one 
energetically open channel (n = 0) and n energetically 
closed channels, where E < hu_\_(2n + 1) for n ^ 0. 
Note that the total amount of the open and the closed 
channels is N and it is related to the principal quantum 
number according to the relation N = n + 1. However, 
the above-mentioned consideration results in an unphys- 
ical asymptotic wavefunction &p(z,p). Equation is 
given in terms of f n (z, p) and g n (z, p), where the motion 
in the z-direction is described by cos q n \z\ and sing n |z| 
(see eq.©). Now, the quantity q n \z\ becomes imagi- 
nary for n (closed channels), and consequently the 
z-dependent terms in eq. ([5j) contain exponentially di- 
verging subterms for z — >• oo. In order to restore the cor- 
rect boundary conditions of the asymptotic wavefunction 
one has to "renormalize " these divergences. This can be 
done within the framework of multichannel quantum de- 
fect theory (MQDT) by partitioning the wavefunction in 
open and closed channel subspaces [33, [3l| . In general 
the wavefunction can be written as follows: 




K 1D K im 

±±-oo ±±oc 
K 1D K 1D 



I, (9) 



where "o" ("c") indicates the open (closed) subspace_re- 
ferring to n = (n / 0), respectively and Y, /, g 
are column vectors dictated by the channel decomposi- 
tion of the considered problem. In order to eliminate the 
diverging terms we choose a linear combination of eq. © 
where Y 00 = I and Y co = -(K 1 ^ - il)' 1 ^. Then eq. 
© takes the following form: 
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where ^ is the "physical" wavefunction which in- 
volves only the open channels since the diverging parts of 
the closed channels are removed. Moreover, the contribu- 
tion of the closed channels during the collision process is 
imprinted in the K-matrix K_l^ phys = K^-\-iK^(I- 



■oo 

I®) provide the 



iKcc^Klo. The roots of det(I - iK 
location of the bound states of the closed channels in the 
parameter space, which energetically lie in the continuum 
of the open channel. 

As it was shown in [Tcj, for collisions which involve a 
single partial wave, the matrix elements of the K 
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K\ D are given by the simple relation {K\ u } 
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(U )n£ t an 5£Ue n >, where after the diagonalization and 



(8) the inversion of (I — i K 
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one can obtain K x ®\ phys 



1D in 



in closed form. However, one should note that K 
this case is a rank one matrix. Interestingly, one can 
rewrite K 1 ®^ phys in a much simpler form by exploiting 

this property with the help of [3 



Then K%> £ phys reads 
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ID, phys 



i 27 



where 7 = Tr(K_l® £ ). Eq. ([IT]) provides us with all the 
transitions which the two atoms undergo during the col- 
lision. Fig. [I] illustrates eq. (fTT|) as a transition diagram, 
where the transitions are not direct, but through the £-th 
partial wave, which participates as an intermediate pro- 
cess. Thus, the first term, both in Fig. [I] and in eq. ([TT]) 
indicates a direct "open-open" (00) channel transition 
and contains, in particular, the description of potential 
resonances. The second term of eq. ([TT]) is related with 
transitions into closed channels, which can be split into 
two processes. As it is shown in Fig. [T]the second term of 
eq. ([IT]) consists of two transition processes, the first de- 
scribes an "open-closed-open" (oco) channel transition, 
while the second describes an "open-closed" (oc) channel 
transition, followed by transitions from all closed chan- 
nels to all closed channels (cc) and finally a "closed-open" 
(co) channel transition. Note that a similar idea has been 
discussed in [33|. 




FIG. 1. A schematic illustration of the transitions into N 
channels which contribute in the scattering process, corre- 
sponding to the if-matrix JC^ phys . The empty (filled) cir- 
cles indicate the open (closed) channels, the arrows depict the 
transitions through the £-th partial wave, and the dots denote 
the N — 3 channels and their corresponding transitions. 

The contribution of the transitions into the closed 
channels during the two-body collision is significantly af- 
fected by the factor 1/(1 — 27), where the parameter 7 
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includes all the closed channel transitions through an ef- 
fective bound state which is supported by all the closed 
channels and for 7 — ^ — i the bound state coincides with 
the threshold yielding a resonant collision process. In- 
deed, according to MQDT this can be justified since the 
parameter 7 = Tr(K_l® £ ) — » — i is a root of the relation 



det(I - iK t 
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1 - iTr(K t 
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and as we have 



mentioned above these roots describe the closed chan- 
nel bound states at given energy E. Thus, we observe 
that the if-matrix approach provides us with an impor- 
tant insight of CIR physics by linking in a transparent 
way the physical picture of virtual excitations Q and the 
Fano-Feshbach-like scheme [3]. Additionally, an impor- 
tant aspect of the transition diagram (see Fig. [T]) is that 
it illustrates in a clear manner the Fano-Feshbach mech- 
anism: the first term represents the continuum (either 
interacting or non-interacting) and the second term ex- 
hibits all the cc transitions, which collectively result to 
an effective closed channel bound state and its coupling 
to the continuum. By applying 7 = — i for each single £ 
partial wave one obtains the condition for resonant scat- 
tering of bosons or spin-polarized fermions under strong 
transverse confinement. Thus for £ = 0,1, 2 one will ob- 
tain the same resonance condition as in 0, S, HH , being 
in agreement with pseudopotential theory. 

Up to now we have been focusing on the case of a single 
partial wave. By extending the if -matrix K_ 3D into an 
£ x £ diagonal matrix, we take into account all the partial 
waves compatible with the atomic species, i.e. even (odd) 
£ for bosononic (fermionic) collisions. Consequently, by 
finding the roots of det(l — iK_l^) of the extended K- 
matrix, we obtain an expression which provides us with 
the positions of all the ^-wave confinement-induced reso- 
nances coupled to each other. 

co oc £—2 

i-ij2Tr[KZ}-- E*«' = °> ( 12 ) 
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£=s +2l'=s 



with = tan ft tan ^ (UuUjv -UjiUw) 2 

00 

and Tr[K£ e }=t!m8 e J2 U ne (13) 
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Note that for eqs. (fT2|) and ([T3]) the summation with 
respect to £ runs over all the even (odd) £ with so = 
(so = 1) for the case of bosons (fermions). Equations ([I2j 
[T3|) represent the main result of the present investigation. 

In order to give some quantitative results in the fol- 
lowing we will focus on low energy bosonic collisions 
(qoa± <C 1) considering s- and d-partial waves. Note 
that a similar analysis, as it is presented below, can be 
done also for the case of spin-polarized fermions which 
involve p- and /- partial waves. 

For the scattering of two bosons the matrix K\® y phys 
can be written in terms of s- and d-wave if-matrices in a 



similar way as in eq. ([TT]) yielding the following relation: 
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where g\ and g^ are defined as follows: 



9i 



1 - iTr\K£ a ] 



(14) 



(15) 



and 
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Tr[KZ}Tr[K^ d ] + Tr[K^ d K^J). 



(16) 



Equation ([T4]) is illustrated by the transition diagram 
in Fig. [2] which depicts all the possible transitions that 
can occur during the collision. In this case there are two 
intermediate transition processes involved, which are re- 
lated to the s- and d-partial wave indicated by solid and 
dash-dotted arrows, respectively. Additionally, one can 
observe in Fig. [2] that the first and second simple terms 
describe direct oo transitions via s- and d-partial waves 
respectively. The third (fourth) term describes cc transi- 
tions via s- (d-) partial waves which are directly coupled 
to the continuum. The last two terms of the diagram 
in Fig. [2] describe processes combining cc transitions of 
both types (s- and d-partial waves-mediated) coupled to 
the continuum. Moreover, in the resonant scattering the 
factors gi and #2 control the contribution of the cc tran- 
sitions via s- and d-partial waves, respectively. One can 
observe that for g\ — » ±oo and #2—^0 the resonant 
scattering is caused by a s-wave CIR, since the cc transi- 
tions through 5-partial wave become the main resonant 
mechanism. Additionally, for #2 —> ioo the resonant 
scattering is caused by a d-wave CIR, since the cc tran- 
sitions through the d-partial wave become dominant, i.e. 
are much larger than the contributions of the transitions 
mediated by the s-partial wave. 

According to eq. ([T2]) the relation which predicts the 
positions of s- and d-wave CIR reads 



i +Cl {k) a - M +om (ad[k) 



+c 3 (fc)- 



m (a d {k) 



0, 



CL±_ \ CL± J CL± 

(17) 

where a s (k) = — tan^ = o(^)/^ 5 a d(k) ~ ~ tanc^=2(&)/& 5 
are the energy dependent s- and d-wave scattering 
lengths and ci(fc), 02(h), cs(k) are constants related to 
the Hurwitz zeta functions [36| , which depend on the to- 
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FIG. 2. A schematic illustration of the transitions into N channels which contribute in the scattering process, for the case of 
s- and d-wave. The empty (filled) circles indicate the open (closed) channels, the solid-black (dashed-dotted) arrows depict the 
transitions through the s- (d-) partial wave, while the dots denote the N — 3 channels and their corresponding transitions. 



tal colliding energy according to the following relations: 
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For go^± <1=> ~ y/2 the values of these constants 
are c x = -1.46035, c 2 = -24.3622 and c 3 = -1.07986. 

Moreover, the "physical" if-matrix, K^?> phys , in the 
open channels reads 
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r>r> 
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(02(A) + C3 (fc)2^1) 
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where the constant C4(fc) is given by the relation: 

c 4 (fc) = 12C(-l/2, 1/2 - (fcaj_/2) 2 ) 

+ (fca ± ) 2 C(l/2,l/2-(fca ± /2) 2 ). 



(20) 



In the limit <2o&± <C 1 the constant C4 takes the value 
c 4 = -5.41533. 

Equation ([T9]) encapsulates all the information of two 
bosons in a harmonic waveguide involving s- and d- 
partial wave scattering. Note that K^' phys possesses 
two singularities and their positions are given by the fol- 
lowing relations: 



a s (k) 



'ci(fc) 



l + ci(fe)^ 



and a d (fe) = * ^— jTT^y fli- 



\ c 2 (/c) + c 3 (fe) 



(21) 

Singularities occur when the 5- and d-wave scattering 
length are approximately half of the length of harmonic 
oscillator a± corresponding to s-wave and d-wave CIR, 
respectively. Furthermore, one observes that the position 
of the d-wave CIR directly depends on the ratio a s (k) / a±_ . 
Interestingly, the second term in the brackets of eq. (p"9]h 
which is related to the d-wave CIR, renders the Fano 
profile of the s- and d- interfering partial waves, where 
the nominator describes the width of the d-wave CIR 
depending strongly on the ratio a s {k)/a±. 

As it was shown in [2[ the scattering of two bosons 
in the presence of a harmonic waveguide can be mapped 
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onto an effective ID Hamiltonian of two bosons which 
interact via an ID delta potential, Vid(z) = giDS(z), 
where giD = h 2 / hclid, with the effective one-dimensional 

2 

scattering length being defined as am = |^( — 1 + 
1.4603^-). The coupling strength gm diverges when 
a>\D —> at the position of the s-wave CIR. This idea 
can be extended in such a way that the d— wave CIR is 



also included in gm giving gm = — — K 
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(22) 



III. RESULTS AND DISCUSSION 

In this section we compare the analytical results of 
the previous subsection with the corresponding numerical 
data. This comparison is based on Cs atom collisions 
representing an ideal system for such a study with a rich 
spectrum of resonances [37], [38| . 

First we evaluate the s- and d-wave energy-dependent 
scattering lengths by solving numerically the free space 
scattering of two Cs atoms interacting via a Lennard- 
Jones potential, V s h(r) = jtt — The dispersion co- 
efficient Cq has been taken from [37j, where the van der 
Waals length is t v dw = (2/iCe/ft 2 ) 1 / 4 = 202 a (a is the 
Bohr radius) and C\i is a free parameter, which controls 
the values of the s- and d-wave scattering lengths. Then 
we use the scattering lengths as an input in eq. ([T9]) , lead- 
ing to a parametrization of the if -matrix K 1 ^' phys in 
terms of C12. The corresponding transmission coefficient 
is calculated by the relation T = 1/[1 + (K_l^ J phys ) 2 ]. 

Moreover, the numerical simulations of Cs-Cs collisions 
in the harmonic waveguide are based on pjl [22|, |39| , 
where we employ the units mc s /2 = H = ojq = 1, with 
mcs being the mass of the Cs atom and ujo = 2ir x 10 
MHz. The longitudinal energy is set to s\\ = 2 x 10 -6 
and the transversal energy is varied within the interval 
2 x 10 -4 < e± < 8 x 10 -4 , corresponding to the ex- 
perimentally accessible range 4tt kHz < uj± < 16tt kHz 
for the waveguide confinement frequency. Additionally, 
the corresponding range of harmonic oscillator length is 
5176 ao < a± < 2588 ao, fulfilling thus the criterion of 
the harmonic oscillator length being much larger than 
the range of the interatomic potential, a±_ ^> £ v dw- 

In Fig. [3] (a,b) we present the analytically and nu- 
merically calculated transmission coefficients, T and T n 
respectively, for several confinement frequencies in the 
vicinity of the d-wave CIR. Notably, the agreement of the 
numerical with the analytical calculations is excellent for 
all confinement frequencies. Moreover, the analytical re- 
sults capture the shift of the <i-wave CIR with increasing 
confinement frequency and the effect of the strong asym- 
metric Fano-lineshape of the transmission spectrum T. 
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FIG. 3. (color online) A comparison plot for (a) analytically 
calculated transmission coefficient T and (b) numerically cal- 
culated transmission coefficient T n of d-wave CIR for several 
confinement frequencies uj±. 



The latter occurs, due to the interference of s- and d-wave 
CIRs, where, as it is indicated in eq. (fT9]h the broad s- 
wave CIR serves as a background and the narrow d-wave 
CIR couples to it. One should note that this difference in 
the widths of the two resonances is attributed to the cen- 
trifugal term, which is absent in the case of s-wave CIR. 
Furthermore, we present in Fig. [4] a high resolution plot 
for two confinement frequencies again comparing the ana- 
lytical with the numerical results. Even in this fine scale 
of the C12 parameter the analytical calculations follow 
the numerical ones both qualitatively and quantitatively. 
The small deviations can be explained by the fact that 
the criterion a± ^> t v dw is not strongly fulfilled for large 
confinement frequencies and consequently the analytical 
results become less accurate. Indeed, one can observe 
that these small deviations in Fig. H](b), which refer to 
the case of strong confinement, become less pronounced 
in the case of intermediate confinement addressed in Fig. 
HI (a). 

Fig. [5] demonstrates the origin of s- and d-wave CIR. 
For that purpose we compare the transmission coefficient 
T (see Figs. [5] (a) and (c)) with the contourplot of the 
expression |det(I— iK\®)\ (see Fig. [5jb)). The transmis- 
sion coefficient T and the relation |det(I — iK\®)\ both 
depend on the total energy E = h 2 k 2 /2fi and the ra- 
tios a s /a_\_, («d/tt±) 5 . More specific, in Fig. [5] (b) for 
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0.999006 0.999008 0.999010 

Cl2 

FIG. 4. (color online) A high resolution comparison of an- 
alytical results (solid line) with numerical data (red squares) 
for the d-wave CIR for two different confinement frequencies, 
uj± = 4 x 10 -4 , 8 x 10 -4 in (a),(b) respectively. 



a given energy equal to the total colliding energy E the 
quantity |det(I — iK\®)\ has two zeros whose position 
is indicated in the contourplot by the dark shaded areas 
(dashed circles), where the positions of the zeros corre- 
spond to two different pairs of ratios (a s /a_\_, (ad/a±) 5 ). 
The latter means that the closed channels possess two 
distinct bound states with energy equal to the total en- 
ergy E at the corresponding values of the pairs (a s /a±, 
(^d/fli) 5 ), where one bound state is related to the s- 
wave ((ad/a±) 5 « 0) and the other one is related to 
the d-wave. On the other hand, in Figs. [5fa) and (c) 
the transmission coefficient T of the open channel at the 
same energy E possesses two zeros which correspond to 
resonant scattering. The position of the <i-wave CIR is 
in Fig. [5fa) at the point a s /a±_ = 0.02 and in Fig. [5fc) 
at the point (ad/a±) 5 = 0.04. Respectively, the position 
of the s-wave CIR is in Fig. [5fa) (see also the inset plot 
on a logarithmic scale) at the point a s /a± = 0.68 and 
in Fig. [5fc) at the point (fld/ai) 5 « 0. Thus, we ob- 
serve an exact matching of the pairs (a s /a_\_, (&d/a±) 5 ) 
at T = 0, denoted by the dotted vertical and horizontal 
lines in Fig. EJb), with the pairs of (a s /a±, (&d/&±) 5 ) 
obtained from the zeros of the relation | det (I — iK\®)\. 
The latter means that the s- and d-wave confinement- 



induced resonances occur due to the corresponding s- 
and d-wave bound states of the closed channels, which 
couple to the continuum of the open channel. Further- 
more, due to the harmonic confinement, s- and d-wave 
CIR are coupled together rendering interference effects, 
i.e. a strong asymmetric Fano-lineshape of the transmis- 
sion coefficient T, as it is depicted in Fig. [5fa) and (c). 
Finally, in Fig. [5fb) the red dashed line indicates the po- 
sition of the (i-wave CIR for any short-range interatomic 
potential according to eq. ([TT]) . 



IV. BRIEF SUMMARY AND CONCLUSIONS 

We have explored two-body scattering within a har- 
monic waveguide taking into account all possible partial 
wave excitations. Our analysis is based on the assump- 
tion that the range of the interatomic interactions has 
to be smaller than the length scale of the harmonic con- 
finement, imposing thus two regimes of distinct symme- 
tries, i.e. spherical for r — >• and cylindrical for r — >• oo. 
The expansion over ^-partial waves leads to an extension 
of CIR physics, where more £ ^ 0-wave CIRs emerge 
being coupled to the broad s-wave CIR predicted in 
Q. Furthermore, we obtain analytically the positions of 
all £-wave confinement-induced resonances (eqs. ([T21IT5|) ). 
This general result refers both to fermionic and bosonic 
collisions. Note that it can be easily extended for more 
than one open channel taking into account inelastic pro- 
cesses as well. Additionally, we have demonstrated that 
the K-matrix approach can provide us with a unique in- 
sight of the underling physics which scattering in confined 
geometries can undergo (eqs. ([TTl dlj) ). More specifically, 
we show that the "physical" ID K-matrix, K^?* phys , 
can be expanded in a Dyson-like form of equations, where 
each term describes all possible transitions which occur 
during the collision process. The latter permits us to 
classify and thoroughly analyze all couplings which are 
induced by the presence of higher ^-partial waves. More- 
over, this scheme allows to unravel the connection of 
the two interpretations for the CIR effect introduced in 
0- El , where all the possible transitions through each im- 
partial wave into the closed channels act in a collective 
way yielding an effective ^-wave bound state supported 
by all closed channels being coupled to the continuum 
of the open channel. The latter yields an asymmetric 
Fano-lineshape in the transmission coefficient T. Addi- 
tionally, we implement the extended analytical model in 
the case of two-body bosonic collisions in the presence 
of a harmonic waveguide with the scattering process in- 
volving only s- and d-partial waves, where we observe an 
excellent agreement between the analytical and numeri- 
cal results. Furthermore, we analyze the d-wave CIR be- 
ing coupled to the s-wave CIR, where both resonances 
depict an asymmetric Fano-lineshape in the transmis- 
sion spectrum. A similar behavior is expected also for 
the case of spin-polarized fermionic collisions involving 
p- and /-partial waves and their couplings. Nevertheless, 
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FIG. 5. (color online) (a) Transmission coefficient T as a function of the a s /a± parameter, (b) contourplot of the quantity 
|det(I — iK\^)\ versus the parameters a s /a± and (ad/a±) 5 and (c) transmission coefficient T as a function of the (ad/a±) 5 . 
In (b) the black dotted lines indicate the position of the resonance with respect to the subfigures (a) and (c), the black dotted 
circles indicate the positions of the zeros of |det(I — iK\®)\ and the red dashed line shows the positions of the resonances for 
an arbitrary short-range interatomic potential. The inset of (a) refers to the dashed boxed area and depicts on a logarithmic 
scale the second zero of the transmission coefficient (s-wave CIR). 



we should mention that for higher partial waves (£ > 4) 
the widths of the corresponding ^-wave CIRs become very 
narrow, due to the pronounced centrifugal barrier which 
the atoms have to tunnel through in order to perform a 
^-partial wave collision. Finally, we would like to mention 
that within the framework of the if-matrix approach the 
origin of the narrow i-wave resonances in the presence 
of confinement is a Fano-Feshbach-like. The latter is not 
in conflict with the interpretation given in [22|, due to 
the fact that when an £ ^ Fano-Feshbach resonance 
is slightly above the threshold, it can be seen also as a 



shape resonance [5(||4l|. This situation occurs also in the 
considered system because the energy of the maximum of 
the centrifugal barrier is three orders of magnitude larger 
than the energy of the ground state of the harmonic con- 
finement, thus the atoms in the <i-wave Fano-Feshbach- 
like state tunnel through the barrier exactly as if they 
would be in a single-channel shape resonant state. This 
complex behavior can be described in a transparent way 
within the if-matrix approach. 



9 



ACKNOWLEDGMENTS financial support by the Deutsche Forschungsgemein- 

schaft. 

The authors are thankful to C.H. Greene and V.S. 
Melezhik for fruitful discussions. P.S. acknowledges the 



[1] L. Tonks, Phys. Rev. 50, 955 (1936); M. Girardeau, J. 
Math. Phys. (N.Y.) 1, 516 (1960); A. Lenard, 7, 1268 
(1966). 

[2] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[3] T. Bergeman, M.G. Moore and M. Olshanii, Phys. Rev. 

Lett. 91, 163201 (2003). 
[4] T. Kinoshita, T. Wegner, and D.S. Weiss, Science 305, 

1125 (2004). 

[5] B. Paredes et al, Nature (London) 429, 277 (2004). 
[6] E. Haller, M. Gustavsson, M.J. Mark, J.G. Danzl, R. 

Hart, G. Pupillo, and H.C. Nagerl, Science 325, 1224 

(2009). 

[7] K. Giinter, T. Stoferle, H. Moritz, M. Kohl, and T. 

Esslinger, Phys. Rev. Lett. 95, 230401 (2005). 
[8] E. Haller, M.J. Mark, R. Hart, J.G. Danzl, L. Re- 

ichsollner, V. Melezhik, P. Schmelcher, and H.C. Nagerl, 

Phys. Rev. Lett. 104, 153203 (2010). 
[9] B. Frohlich et al, Phys. Rev. Lett. 106, 105301 (2011). 
[10] B.E. Granger, and D. Blume, Phys. Rev. Lett. 92, 133202 

(2004). 

[11] J.I. Kim, V.S. Melezhik, and P. Schmelcher, Phys. Rev. 

Lett. 97, 193203 (2006); Rep. Progr. Theor. Phys. Supp. 

166, 159 (2007). 
[12] V.S. Melezhik, and P. Schmelcher, New J. Phys. 11, 

073031 (2009). 

[13] D.S. Petrov and G.V. Shlyapnikov, Phys. Rev. A 64, 

012706 (2001). 
[14] L. Pricoupenko, Phys. Rev. Lett. 96, 050401 (2006). 
[15] Z. Idziaszek and T. Calarco, Phys. Rev. Lett. 96, 013201 

(2006). 

[16] V. Peano, M. Thorwart, C. Mora, and R. Egger, New J. 

Phys. 7, 192 (2005). 
[17] S. Sala, P.-I. Schneider and A. Saenz, Phys. Rev. Lett. 

109, 073201 (2012). 
[18] SG. Peng, H. Hu, X.J. Liu and PD. Drummond, Phys. 

Rev. A 84, 043619 (2011). 



[19] S. Saeidian, V.S. Melezhik and P. Schmelcher, Phys. Rev. 

A 77, 042721 (2008). 
[20] V.S. Melezhik and P. Schmelcher, Phys. Rev. A 84, 

042712 (2011). 

[21] Y. Nishida, and S. Tan, Phys. Rev. A 82, 062713 (2010). 

[22] P. Giannakeas, V.S. Melezhik and P. Schmelcher, Phys. 
Rev. A 84, 023618 (2011). 

[23] A.M. Rey, R. Sensarma, S. Foiling, M. Greiner, E. Dern- 
ier, and M.D. Lukin, Europhys. Lett. 87, 60001 (2009). 

[24] B. Deb and J. Hazra, Phys. Rev. Lett. 103, 023201 
(2009). 

[25] W. Hofstetter, J.I. Cirac, P. Zoller, E. Dernier, and M.D. 

Lukin, Phys. Rev. Lett. 89, 220407 (2002). 
[26] B. Deb. larXiv: 1205.0644^ 2. 

[27] DA. Harmin, Phys. Rev. Lett. 49, 128 (1982); Phys. 

Rev. A 26, 2656 (1982); Comments At. Mol. Phys. 15, 

281 (1985). 
[28] U. Fano, Phys. Rev. A 24, 619 (1981). 
[29] C.H. Greene, Phys. Rev. A 36, 4236 (1987). 
[30] M. Aymar, C.H. Greene, and E. Luc-Koenig, Rev. Mod. 

Phys. 91, 1015 (1996). 
[31] PG. Burke, R-Matrix Theory of Atomic Collisions, 

Springer- Verlag, Berlin, Heidelberg, 2011. 
[32] K.S. Miller, Mathematics Magazine 54, 67 (1981). 
[33] G. Alber and P. Zoller, Phys. Rep. 199, 231 (1991). 
[34] A. Imambekov, A. A. Lukyanov, L.I. Glazman, and V. 

Gritsev, Phys. Rev. Lett. 104, 040402 (2010). 
[35] P. Giannakeas, V.S. Melezhik, and P. Schmelcher, Phys. 

Rev. A 85, 042703 (2012). 
[36] Handbook of Mathematical Functions, M. Abramowitz 

and I.E. Stegun (Department of Commerce, Washington, 

DC, 1972), 10th ed. 
[37] C. Chin et al, Phys. Rev. A 70, 032701 (2004). 
[38] M. Mark et al, Phys. Rev. A 76, 042514 (2007). 
[39] V.S. Melezhik, Multi-Channel Computations in Low- 
Dimensional Few-Body Physics, Springer- Verlag, Berlin, 

Heidelberg, 2012. 
[40] B. Gao, Phys. Rev. A 80, 012702 (2009). 
[41] B. Gao, Phys. Rev. A 84, 022706 (2011). 



